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Abstract. We investigate theoretically the combination of first-order quadrupole-quadrupole and second- 
order dipole-dipole effects on the long-range electrostatic interactions between a ground-state homonuclear 
alkali-metal dimer and an excited alkali-metal atom. As the electrostatic energy is comparable to the dimer 
rotational structure, we develop a general description of the long-range interactions in the framework of 
the second-order degenerate perturbation theory, which allows for couplings between the dimer rotational 
levels. The resulting adiabatic potential energy curves exhibit avoided crossings, and cannot be expanded 
on the usual 1/R n series. We study in details the breakdown of this approximation in the particular case 
Cs 2 +Cs(6 2 P). Our results are found promising to achieve photoassociation of ultracold trimers. 
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1 Introduction of both external and internal degrees of freedom of the 

molecules offers amazing opportunities for fundamental 
Researches on cold (T < 1 K) and ultracold (T < 1 mK) and applied research, like high-resolution molecular spec- 
gases have attracted a considerable interest over the last troscopy [3, 4, ; 5.„.&] , resonant dynamics [THS] , phase-driven 
years, as at such temperatures, the properties of the gases chemistry [TDJ[TT], quantum degeneracy [T^HSIE] , quan- 
are predominantly governed by purely quantum effects, turn phase transitions [THIUSE] ' universal quantum states 
like the tunneling, quantum resonances, and quantum Bose dEEMO] , metrology of fundamental constants [2TJ1221I231 
and Fermi statistics. One of the most recent developments 24J, quantum information [25l[26l[27] , cold and ultracold 
concerns cold and ultracold molecules - the central theme collisions and chemistry p51EM3l l3TTl3^1 [3^1[5il[33] . 
of the present topical issue [UH]. The exquisite control 



2 M. Lepers, O. Dulieu: Cold atom-molecule photoassociation: long-range interactions beyond the 1/7?" expansion 

The central feature of cold and ultracold gases is that the BC molecule. In contrast with the atom-atom case, 
the relative kinetic energy of the particles is so tiny that the long-range interactions between atoms and molecules 
their interactions is dominated by long-range forces in- has been less often addressed quantitatively in the liter- 
duced by the spatial extension of their charge distribution, ature. In most studies, they have been evaluated at fixed 
The related potential energy is usually described by the geometries in order to match or to fit at large distances 
well-known multipolar expansion ^ n C n / R n in inverse potential energy surfaces obtained by quantum chemistry 
powers of the (large) distance R between the interacting computations [60,61,62,63,64]. Note however that in a re- 
particles [36,37,38]. With the development of ultracold cent paper [55] > Kotochigova reported a study on the van 
physics, a vast amount of accurate calculations of the C n der Waals interaction between a ground state alkali metal 
coefficients is available nowadays for several atomic species atom and a ground state alkali-metal diatomic in a given 
like alkali-metal [39,40,41,42,43,44], alkaline-earth [45], rovibrational level, i.e. beyond the fixed geometry frame- 
rare gas atoms [35] , ° r open-shell atoms |4"71l48] . work. As an illustration, we focus here on the PA reaction 



In the present study, as the third of a series of papers 
|49ll50] (hereafter referred to as papers I and II, respec- 
tively) , we investigate the long-range interaction of an ex- 
cited alkali-metal atom and a ground state alkali-metal 



between an ultracold ground state Cs2 molecule in a vi- 
brational level Vd and a rotational level N and a cesium 
atom 

Qs 2 {X 1 S+ v d , N) + Cs(6 2 5) + hv -> Cs^ (1) 



diatomic molecule. Indeed, samples of ultracold ground which could be explored in the experimental group in Or- 
state alkali-metal diatomics are available in several exper- say. In paper I we have determined the leading term of the 
iments, surrounded in most cases by ultracold atoms [511 long-range interaction, namely the quadrupole-quadrupole 
I52TI55 , 54,55,56,57. 58] . Tuning the frequency v of a laser term between an excited Cs(6 2 P) atom and a ground 
to the red of a resonant atomic transition, it should be state Cs2(A 1 Z'+, Vd, N) molecule, with N > 0. We ex- 
possible to induce the photoassociation (PA) of an atom tended this study in paper II to the calculation of van der 
A and a molecule BC according to A + BC + hv — > (ABC)* Waals coefficients with second-order perturbation theory, 
just like it is routinely achieved for atoms [59] . The excited Due to the competition between the rotational energy of 
complex ABC* could spontaneously decay afterwards in a the molecule and the electrostatic interactions, we demon- 
stable ultracold trimer ABC, providing then an example strated that the validity of the multipolar expansion at 
of a photo-assisted ultracold chemical reaction. While the short distances is not limited by the overlap between elec- 
latter process most probably relies on complex chemical tronic clouds of the particles (occurring around 45 a.u. 
forces, the former PA step is obviously controlled by the in the present case), but by the degeneracy between vari- 
long-range interaction between the excited atom A* and ous long-range states of the complex occurring in the 80- 



M. Lepers, O. Dulieu: Cold atom-molecule photoassociation: long-range interactions beyond the 1/R n expansion 



100 a.u. distance range (1 a. u.=ao=0. 0529177 nm). The 
aim of this article is therefore to investigate the distance 
range for which the multipolar expansion is still applica- 
ble to compute long-range interactions, but resulting into 
potential curves determined by a second-order degener- 
ate perturbation theory which cannot be written anymore 
as the usual ^2 n C n / R n expression. The paper is orga- 
nized as follows. In Section [21 we recall the basics of the 
perturbative calculations of long-range interactions, as de- 
scribed in papers I and II, and reformulate it up to the 
level of the second-order degenerate perturbation theory. 
In Section [21 we derive the expressions of the first-order 
and second-order interaction energy. The results for the 
interaction of a ground state Cs2 molecule and an excited 
Cs atom are presented in Section [H as the initial state of 
the photoassociation of a cold atom-molecule pair, before 
addressing some prospects in Section [51 

2 Energy scales and perturbation theory 

In our two previous articles I and II, the long-range in- 
teractions between a homonuclear alkali-metal dimer and 
an excited alkali-metal atom is investigated, choosing the 
example of Cs2+Cs*. The dimer is in its ground electronic 




Fig. 1. The two coordinate systems, XaYaZa (D-CS) and 
XYZ (T-CS) defined for the dimer and for the trimer, respec- 
tively. The Za axis is along the dimer axis, while Z is oriented 
from the center of mass of the dimer towards the atom B. The 
Y and Ya axes coincide and point into the plane of the figure. 
The subsystem A in this figure is the Cs2 molecule, the subsys- 
tem B is the Cs atom. The T-CS is related to the laboratory 
coordinate system (xyz) by the usual Euler angles (a,/3,"/), 
not represented here. 

axis. The fine structure of the atomic level will be included 
in a further study. 

We started from the general form of the electrostatic 
energy between two interacting charge distributions A (the 



state \X~>-E+), in an arbitrary vibrational state |u d ), and dimer ) and B ( the at0 m), separated by the distance R, 
rotational state characterized by the quantum number N wr itten as the well-known multipolar expansion: 



with projection m on the Z axis of the trimer coordinate 
system (T-CS) (Fig. [TJ). The alkali- metal atom is in its 
first excited state 2 P characterized with general indexes 



1 



v el { R )= J2 E 

L A ,L B =n M=-L < 



ftl+L A +L B 



(2) 



labeling the valence electron: n for the principal quantum where L < is the minimum of La and Lb- Each multipole 
number, I for the orbital one and its projection A on the Z of order Lx (with X — A, B) is associated to the tensor 
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operator QfJ x (fx), which can be expressed in a coordinate ting La = Lb = 2 in Eq. Q 
system whose origin is the center of mass of X: „ 24 ^ Qif (r^)Q^ M (fs) 



'"' i? 5 4^ (1 + Af)l (1 - Af)! ' 1 ' 



Ql'x (fx) — J -zj — — J]] <kff x Ylx (@i> ' (3) The corresponding potential energy curves BqN(N + 1) 



where is the value of each charge i composing X, and 
Yjf x are the usual spherical harmonics. In Eq. ([2]), the 
assumption has been made that the quantization axis is 
Z, oriented from the center of mass of A to the center of 
mass of B, hence the factor /l a l b m reads 



C$1 i? 5 as functions of R have been reported in paper I. In 
the case of Cs2+Cs*, the curves start to cross each other 
at distances R m « 100 a.u.. This distance is higher than 
the Leroy radius [07] estimated at Rlr ~ 45 a.u. where 
electronic clouds start to overlap. Therefore we define the 
crossing region by Rlr < R ^ R m where the perturbation 



i \\^ B ( t 4-7, V 

f T T = I ! \ A ~r Bh approach described previously is not applicable. This is a 

TLaLbM y/(L A + M)\(L A -M)[ 

1 general feature due to the competition of the rotational 

(4) 



\J(Lb + M)- {Lb M). energy of the dimer with the electrosatic energy. Note that 

In Papers I and II, we defined the zeroth-order energy £ ® for Li 2 +Li, we estimate this crossing in the 26 < R < 
of a given (2N + 1)(2£ + 1) times degenerate state labelled 43 a.u. range. 

p as In paper II, we studied the second-order contribution 

c o i tt\ /ri to the dipole-dipole interaction, i.e. La = Ln = 1. scaling 

as R~ 6 . We calculated the Cq coefficient corresponding to 
each eigenvector \<fp) of V qq in the subspace of degeneracy 
associated to £° Noting that 



where E n ( is the energy of the free atom, Ex Vd N = Ex Vd + 
Bxv d N (N + 1) the free energy of the dimer in the rovi- 
brational level (vd, N) considered for simplicity sake as a 
rigid rotor with a rotational constant Bxv d (Bxv d =o = 
B = 1.173 x 10~ 2 cm -1 [66J). In the following, the en- we write the Cq coefficient in the general form 
ergy origin will be set at Ex. Vd =o + E^p — without loss 
of generality of our treatment. As the two reactants get 



7} I A 



C 6 = R 6 ^2 X! C rni\iCm 2 \ 2 (™lAl| Vj/ |m 2 A 2 ) , (8) 
miAi rn 2 Ao 



closer from each other, the degeneracy over the different where ^dd is the second-order dipole-dipole operator 

sublevels is progressively lifted, and the electrostatic en- y( 2 ) _ _4 \ 

M ^w' { > Ea ~ Exv * N) + { - Ei ~ Ent) 

ergy must be calculated using the perturbation theory for 

x gyjgo) ($a\ Qi M Qi M \$b) W2f ( ) 

degenerate levels. (1 + M)\ (1 - M)\ (1 + M')\ (1 - M')\ ' 1 ' 

The first-order correction is the quadrupole-quadrupole where a = X'v' d N'm' and b = n'l'X 1 are the quantum 
energy of interaction V qq scaling as R~ 5 obtained by set- numbers associated to the intermediate states accessible 
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by dipolar transition. In paper II, the calculations for pect those interactions to be significantly smaller (at least 
Cs2+Cs* revealed that, in the crossing region, the second- one order of magnitude), in the range of distance that we 
order dipolar energy becomes comparable to the first-order consider here, 
quadrupolar and the rotational energies. 

In order to describe the electrostatic interactions in 3 Calculation of the electrostatic energy 
the crossing region, we need to formulate the perturba- 

This section is devoted to the calculations of the matrix 

tion approach in a different way. As the rotational and 

aI ATVl/Ml + in S~\ T I 1 1 1 1 1 T 



elements of V qq and V dd in the basis {|7VmA)}. The prin- 



electrostatic energies are comparable, we consider the per- 

ciple of the calculations has been presented in Papers I 

turbation hamiltonian W 

and II. So we will only recall the main steps, emphasizing 



.2 ... M 2) 



W = B Xvd N +V qq + V£>. (10) 

HH aa on the most important difference: the couplings between 

so that the zeroth-order energy of a state p is reduced to distinct values of the rotational quantum number N are 

Ep = E Xvd + E n g. Each value of E® is associated to a taken into account. 

subspace of quasi-degeneracy whose basis vectors can be 

labeled by N, m and A quantum numbers. In the present 3.1 First-order quadrupole-quadrupole interaction 

formulation, the electrostatic interactions now couple dif- 
ferent values of N, so that each subspace of degeneracy 
has in principle an infinite dimension. 



("2") 

It may seem surprising that, in Eq. (|10p. V dd is actu- 
ally manipulated like a first-order contribution, whereas it 
is a second-order term. To justify Eq. (|10p. we have cal- 



As the atomic part of the matrix elements of V qq are un- 
changed in comparison to Eq. (13) of Paper I, we only 
focus here on the dimer part. The wave function !f r Ar m (5) 
associated to the rotational state \Nm) of the dimer is 
written in the T-CS as 



+ 1 , N 



cola.ed the second-order correction to energies and state ^ = V —*•«>. <"> 

vectors for degenerate unperturbed levels, and found that where d^ l0 (5) is the reduced Wigner matrix element. The 



the corrections to energy are given by the eigenvalues prefactor yj (27V + l)/2 ensures normalization to unity with 
of the operator V^f . In principle, in order to calculate respect to 5. Therefore, we can write the matrix element 
the total second-order contribution of the electrostatic en- of the quadrupole moment of the dimer as 



ergy, dipole-quadrupole, quadrupole-dipole, quadrupole- I N m 



N 2 m 2 



quadrupole, ... terms must be taken into account. However a/(2/Vi + 1)(2A?2 + 1) 

= 2 9 ° 

one should keep in mind that those interactions scale as 



x / dS S mSdZ\ (S)d 2 MO (S)dZl (S) 

R , R and R respectively. If we refer to the values 



of C$ and C10 coefficients (see e.g. Ref.[40j), we can ex- 



/ 2A 2 + l ^ATtO r N imi (-19 s ! 
V 2iVi + 1 N ^20^ N 2 m 2 2M ( i2 ' 
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where q 2 is the (only non-zero) tensor component of the where a denotes in a short way the intermediate state 
quadrupole moment of the dimer in its coordinate system \X'A'v' d N'm'X'). For the sake of clarity, we explicitly in- 
(the D-CS). In Eq. (fT2)) . we have written the quadrupole troduce A and yl',the projection of the dimer orbital mo- 
moment operator (see Eq.© with L — 2) using the prop- mentum on the Za axis (A — in our particular case). The 





ertyU e~^d L M0 (9) = ^4tt/(2L + l)Yj 
have calculated the integral of three Wigner functions 



Then we quantity x{a) is part of the dimer polarizability ol^I^? 



and obtained the Clebsch-Gordan coefficients C^ b g. So 
the matrix element of V qq reads 



(N x mxXi V qq N 2 m 2 X 2 ^ 



4 (rl t ) 



- , AT 2 020^^020 
Lj N 2 m 2 2M [ ~ y l\22-M 

M=-2 



(2 + M)\{2- M)\ 



(13) 



We recall from Paper I, that for Cs2, we estimated q® 
at 18.58 a.u., and we calculated (f*g P ) = 62.65 a.u. for 
cesium using an Hatree-Fock method. Equation (|13[) im- 
pose strong selection rules: (i) the projection of the total 
orbital quantum number mj = m + A is conserved; (ii) 
N\ = N2,N2 ± 2, which means that the parity of N is 
conserved. 

3.2 The second-order dipolar interaction 

Once again, only the contribution from the dimer is modi- 
fied with respect to the developments of Paper II. In order 



to calculate (NxmiX 



V, 



(2) 



,1,1 



N 2 m 2 X 2 /, we concentrate on 



the calculation of the quantity 



x(a) = (XAvdNiiniAi 



X'A'v' d N'm'\' 



x (X'A'v' d N'm'X' 



Q^ M ' XAv d N 2 m 2 X 2 ) , (14) 



1 In our case, 8 — 5, whereas (f> is not defined so that we take 
i = 0. 



(see Eq. (12) of Paper II) according to the relation 



CffW = 2(-D M E £r (15) 



where z is a real or imaginary frequency. The dipolar ma- 
trix element is 



XAv d Nm 



= ]T (Nm\d^\N'm l ){XAv d \%\X'A'v' d ), (16) 

where fi represents the different possible projections on 
the Za axis, and the corresponding tensor components 
of the dimer dipole moment, fi = is associated to £ — > £ 
transitions and [i = ±1 to £ — > II transitions. In a similar 
way as Eq. (fl"2"|) , and by use of the property [SB] 



CZbp = (-1)° 



we write the rotational part 



2c+l^ b(3 



26+1 



a 



cya— ct' 



(17) 



{Nm\d\^\N'm') 



( i\M+/j J 2N +1 rN ' m > r N'-ii n %\ 

V 2N' + 1 Nml ~ M noi-h- 



Finally, Eq. (IM|) becomes 



_ M V /(2jV 1 + l)(2iV2 + l) 

-w-i L ) 2N , + 1 

/-iN'm' /-iN'm' \ y /-tN' — fi 

x iV 1 mil-M iV2m2l-M' / ^NiOl-u 

* C»'-f_ ll \{XAv d \q>£\X'A'v d )\\ (19) 



and 
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polarizability, [i n 'f,ni = r n 't>, neC^ Q ° w /V3 is the atomic 



'MM 



AT' 
X C 



V(2iVi + 1) {2N 2 + 1) 



2N' + 1 

N'm' 

N m /~iN'm 
Nxmil-M^ N 2 m 2 l-M' 



o iVi010 o JV2010 a || V Z J 
+ 2Cwi011 C/VaOll^-L (z) 



(20) 



transition dipole moment, AE n i^^ n i — E n >^ — E n g is the 
atomic excitation energy and Q(x) is Heaviside's function, 
here equal to 1 if AE n ig^ n g < 0. 

Equation (|2ip indicates to us that the selection rules 

" (2) 

associated to V dd are the same as V qq . Therefore the sub- 



with a || and a± respectively the parallel and perpendic- s P ace of degeneracy associated to a given zeroth-order en- 

ular dipole polarizabilities of the dimer in the X elec- K can be divided illto independent subspaces de- 

tronic state and v d vibrational level (see Paper II for a fined: « b ^ the Projection of the total orbital momentum 

description of the polarizability calculation for Cs 2 and ™j = m + A > and (") h Y the P arit y of the rotational quan- 



Cs). Therefore, the matrix element of V dd looks pretty 
much like Eq. (32) of Paper II: 

Nim 1 X 1 \v d i f\N 2 m 2 X^ 

^ ^ ^ (l + M)\ (l - MY. (1 + M'Y. (1 - M'Y. 

MM' N'm' i'\' v ' y ' y ' y ' 



^{2N X + 1) (2N 2 + 1) , m . 

'^Nm, 1 — M^l\ 



21- 



E 



2N 

A' s~i 

u 'AilM°' 



l 



1 -rf'A' s~<£'\' 



2£/ _|_ " y fAilM v - y £A 2 lM' 



-"TVmil-M^JVmal-M' 



-»iV'0 



ii 



^ Ly JVi010 L/ Af 2 010 t 



+2C^'J )n C^'J )n a±(iu}) I otntn't'^tu) 



-40{-AE n 



riN'o 

TVi010 ( - / Ar 2 010 



(AE n 



a\\ (AE n w,ni) 

2 



(m 



n'l'nt) 



f EE - 

M AT'm' 



1 



[(1 + M)!(1 -M)!] 



7(2^ + 1) (2JV 2 + 1) C * 



2iV 
duj 



Nm 1 lM^Nm 2 lM in 1 m2 \i\2 



N'i 



( f-iN'O 
I °JViO 



OlO^JVaOlO"!! 



JV'O 



an [iw) 



(21) 



turn number N. 

4 Results and discussions 

The first-order correction to energy due to the operator 
W is obtained after a diagonalization in each subspace 
of quasi-degeneracy associated to E® for each value of R. 
The eigenvalues E* (R) define potential energy curves with 
avoided crossings, that will thus be called adiabatic. The 
corresponding eigenvectors will be labeled |(p) \m,j \ ,N), 
where \mj\ = |to + A|, N is the dimer rotational quan- 
tum number to which the state tends as R — > oo, and p 
a number starting from 1 and increasing with energy in a 
given symmetry (that is to say a given \mj\ and parity of 
N). By contrast, the potential energy curves and eigen- 
vectors obtained in Papers I and II are called diabatic, as 
they cross inside in a given symmetry. The eigenvectors 
are labelled (p) \mj\ , N d ), the d superscript standing for 



In this equation a n t n 't' is the state-to-state polarizability diabatic. Unlike the adiabatic states, the diabatic ones are 
corresponding to the n£ — > n'£' atomic transition (see Eq. totally included in the manifold defined by a given value 
(22) of Paper II), a c is the core contribution to the atomic of the dimer rotational quantum number N. 
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Fig. 2. Long-range potential energy curves between a ground- 
state Cs2 and an excited Cs(6 2 P) atom, as functions of their 
mutual separation R, for the S + symmetry and for: (a) the 
even values of JV and (b) the odd values of N. The range 
of energy corresponds to the rotational energy of the dimer 
B N(N + 1) for JV = to 5. The adiabatic curves, result- 
ing from the diagonalization of W (see Eq. (|10p ). are plot- 
ted in full lines. The diabatic ones, given by BoN(N + 1) + 
C a /R 5 + C 6 /R 6 (see Paper II) are plotted in dashed lines. The 
curves in heavier lines are those involved in the crossing be- 
tween \(3)E + ,N = 2 d ) and \{4)E + ,N = 4 d ), discussed in the 
text. 

Both the diabatic and adiabatic potential energy curves 
are shown on Figs. [2HH in dashed and solid lines respec- 
tively. Each panel corresponds to given \mj\ and parity 
of N. Following our previous work, we have labelled our 
energy curves with the diatomic-like symmetries S , II, 
A, <!>, r and H. All the curves displayed on Figs. ©-© 
present similar general features. On their right part, for 
R > 100 a.u., the diabatic and adiabatic curves are not 
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Fig. 3. Same as Figure [5] except that the range of energy 
corresponds to the dissociation limits characterized by N — 8 
to 12. The curves in heavier lines are those involved in the 
crossing between |(ll)X' + ,iV = 10 d ) and |(12)X' + ,iV = 12 d ), 
discussed in the text. 
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Fig. 4. Same as Figure [5] but for the S~ symmetry. 

distinguishable. This corresponds to the region of space in 
which the usual perturbation approach, characterized by 
the inequalities 



B Xvd N 2 ) » (V qq ) » (V$) , 



(22) 



is applicable. In the central part, for 50 < R < 100 a.u., 
strong differences are visible between the diabatic and adi- 
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Fig. 5. Same as Figure [2] but for the 77 symmetry. 
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Fig. 6. Same as Figure [2] but for the A symmetry. 
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Fig. 7. Same as Figure [2] but for the $ symmetry. 
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Fig. 8. Same as Figure [2] but for the T symmetry. 



abatic curves. In particular, we can see avoided crossings ues of ^ each subspace of qua si-degeneracy has strictly 

of the adiabatic curves, whereas the diabatic ones do cross. speaking an infinite dimension. Of course, in our calcula- 

At last, in the left part, for R < 50 a.u., we see a dense tionSj we haye chosen a maximum value of N ^ /v mox) and 

set of attractive curves connected to higher-TV dissociation haye tested thg convergence with respect to that value . 



limits. 



Our test indicate that, in order to describe properly the 



Since the quadrupolar and dipolar parts of the Hamil- energy curves up to an energy of BqN*(N* + 1), one has 



tonian W (see Eqs. (|13l) and (J^TJ) couple different val- to chose at least N max N* + 4. For our calculations, 
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Fig. 9. Same as Figure [5] but for the H symmetry. 

we have taken N max = 17, which gives us a reliable de- 
scription of the range of energy in Fig. [3] This observation 
can be explained as follows. In the range of distances in 
which we are interested, there exists a value of R where 
the most attractive curve coming from the TV* + 2 dissoci- 
ation limit has the energy R>oN*(N* + 1) corresponding to 
the TV* dissociation limit. So, in order to correctly describe 
that TV* + 2 curve, one has to include in the calculations 
the curves connected to the TV* + 4 dissociation limit, as 
they are mutually coupled. By ignoring the Ce/R 6 part 
of the energy of the TV* + 2 curve, we can estimate that 
the crossing region scales as R ~ (42V + 6) 1/5 . This weak 
sensitivity to TV explains why the difference N max — TV* 
does not vary significantly with increasing TV*, although 
the rotational levels of the dimer get further away from 
each other. This is particularly the case for the low ener- 
gies considered here, which are characterized by TV* < 20. 

It is worthwhile noting that the previous reasoning 
does not hold for R < 40 a.u.. In this region, the con- 
vergence on N max is very difficult to obtain. Due to the 
very strong independence of the electrostatic interaction, 



the influence of the TV* ± 2 curves is not the only one 
significant on the TV* curves. We have not represented 
such a region on Figs. [2HS1 since the electronic overlap 
will certainly play an important role. However, it cannot 
be excluded that such a convergence problem arises with 
different reactants. 

In the region 40 < R < 100 a.u., the interactions be- 
tween different diabatic curves strongly modify the as- 
pect of the potential energy curves. The most obvious 
modifications are the clear avoided crossings of the adi- 
abatic curves. This is for example the case between the 
\(3)S + ,N = 2 d ) and the \(4)E+,N = 4 rf ) states, or be- 
tween the | TV = 10 d ) and the | (12) S+, TV = 12 d ) 
(see Figs. ([2]) and ((3]) respectively). More generally, for 
the S + symmetry, there is an avoided crossing between 
the less attractive curve of the TV manifold and the more 
attractive one in the TV + 2 manifold. In non-17 symme- 
tries, the higher density of states (3 for each TV instead of 
2 for S + ) alters the usual two-state picture of the Landau- 
Zener problem, due to the influence of neighboring curves. 
At last, for the £~ symmetry, we cannot see any Landau- 
Zener crossing, as the diabatic curves are less numerous 
and approximately parallel. 

We have estimated the Landau-Zener probability P pr 
to go through the crossing between the curves (p) and 
(?') in a diabatic way, and applied the result to the two 
crossings evoked above. We start from the general formula 



P pr = exp (~2nr pr ) , 



(23) 
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with 

r, 



pr v p (Ro,T) \^W* p (R = R ) - ^W r d r (R = Ro)\ ' 

(24) 

In Eq. (|2l?|) , W pr (R) is the matrix element of the hamilto- 
nian (jlOp in the diabatic representation 

W* r (R) = (ip)mnN<t \w(R)\ (r)m. J2 N*) , (25) 

Ro is the distance at which the crossing occurs, and v(R, T) 
is the temperature-dependent classical velocity of a parti- 
cle entering the crossing in the p channel, 



v P (R,T) 



! 2(W PP (R^ oo) - W pp (R)) + 3k B T 
A 4 



(26) 



with k B Boltzmann's constant and fi the atom-dimer re- 
duced mass. The off-diagonal term W pr (R) appearing in 
Eq. (f2"3")l can be written 



u 5 i °6 



(27) 



where C 5 and C 6 are the off-diagonal matrix elements: 



C 5 = R 5 (p 



V, 



QQ 



and C 6 = R 6 (p 



V, 



(2) 



<I<1 



derivative -§rW pp (R) reads 



5C 5 6C 6 
i? 6 R7 ' 



r); and the 



(28) 



entrance channel p 


v p (10 _6 xa.u.) 


Ppr (96) 


\(3)E+,N = 2 d ) 


2.4 


4 


(4)r+,iV = 4 d ) 


4.9 


21 



Table 1. Characterization of the Landau-Zener crossing 
between the , N = 2 d ) and | (4)X" + , TV = 4 d ) diabatic 

curves, with respect to each entrance channel. The crossing 
occurs at Ro = 59.3 a.u., and the corresponding coupling W pr 
is 0.058 cm -1 . The velocities v p are calculated for a tempera- 
ture T = 1 mK. 



entrance channel p 


v p (10" 6 xa.u.) 


Ppr (%) 


\(U)E+,N = 10 d ) 


3.6 


89 


\(12)S+,N = 12 d ) 


8.6 


95 



We have applied Eq. ([2"5]) to the crossings between 
\(3)£ + ,N = 2 d ) and \(4)E+,N = 4 d ), and between | (11)^7 
and |(12)Z ,+ , N — 12 d ), and we have summarized the re- 
sults in Tables Q] and [2j respectively. The dynamics can 
drastically change from one case to another, turning from 
almost adiabatic in the lower channel of the former cross- 
ing (P pr = 4%), to almost diabatic in the latter crossing 
(Ppr —95%). Obviously, for a given crossing, the dynam- 
ics of a particle coming from the upper channel is more 



Table 2. Characterization of the Landau-Zener crossing be- 
tween the \(11)E+,N= W d ) and \(12)S + ,N = 12 d ) diabatic 
curves, with respect to each entrance channel. The crossing oc- 
curs at Ro — 47.6 a.u., and the corresponding coupling W pr is 
-0.027 cm -1 . The velocities v p are calculated for a temperature 
T = 1 mK. 

diabatic, as it has acquired more kinetic energy at the 
crossing point. However, in the crossing between the (11) 
and (12) curves, even if the velocity of the upper channel is 
more than twice the velocity of the lower one, the Landau- 
^ffliez Jfls^ability only increases by 6%. Actually, the most 
determinant parameter for the nature of the crossing is the 
coupling matrix element W pr . The second crossing tends 
to be more diabatic, because | W pr | is significantly smaller, 
i.e. -0.027 cm -1 with respect to 0.058 cm -1 . This shrink- 
ing is due to the competition between C 5 > and C 6 < 0. 

Apart from the Landau-Zener-like crossings, the dia- 
batic curves can also change their behavior in a "smoother" 



12 M. Lepers, O. Dulieu: Cold atom-molecule photoassociation: long-range interactions beyond the 1/7?" expansion 



way as a function of R. This is for example the case for 
the and (1)77 states correlated to N — 0, and 

for the (2)S + and (2)77 states correlated to N = 1. For 
7? > 100 a.u., those states are characterized by C5 = 0; but 
as 7? decreases, they progressively acquire a R~ 5 behav- 
ior, due to the coupling with the higher curves. This result 
might be of strong importance, especially in the prospect 
of creating an excited trimer from a N = sample of 
dimers. 

This variety of coupling situations has made us search- 
ing for a better criterion of applicability of the C n /R n ex- 
pansion, than the one proposed in Paper I, and based of 
the curve crossings. Basically, we estimate that the C n /R n 
expansion stops to be valid as soon as the ratio between 
the coupling and the energy difference between two states 
exceeds a certain threshold e. Therefore we define a state- 
to-state normalized coupling W pr as 



_ W* r (R) 
Wpr{R) ~ W*{R)-W*.(R)' 



(29) 



where Wp r (R) are the matrix elements of W in the di- 
abatic representation, and we consider that the diabatic 
representation fails for 7? such that 



max I Wpj-fiZ)] > e. 



(30) 



In Eq. (|30p . the equality corresponds to the threshold dis- 
tance R = 7?*, which depends on the state \{p)mjN d ') 
under consideration. The quantity W pr is plotted on 
Figs. [Tt)lll2l for various symmetries. In all cases, the quan- 
tities W pr tend to zero with increasing atom-dimer dis- 
tances (note that 7? goes up to 500 a.u. on Fig. ITDlfTS|) . 
since the off-diagonal elements of W vanish. On the con- 



1= 



(2)E + ,(V=2-(1)i: + ,iV=0 
(2)E + ,jV=2 - (3)E + ,iV=2 
(2)Z + ,/V=2 - (4)E + ,iV=4 
(2)E + ,jV=2 - (5)i: + ,iV=4 




500 



(3)E + ,/V=2-(1)5: + ,A<=0 
(3)Z + ,/V=2 - (2)S + ,iV=2 
(3)S + ,(V=2 - (4)E + ,iV=4 
(3)£ + ,/V=2 - (5)E + ,iV=4 




500 



Atom-dimer distance R (a.u.) 



Fig. 10. Normalized state-to-state couplings W pr (see Eq. 
(|29|l 'l as functions of R, for the two diabatic states of the 
£ + ,N = 2 manifold: (a) state (2), (b) state (3). The two hor- 
izontal dotted lines correspond to IWprl = e — 0.1. 



(2)n,jv=2 - (i)n,jv=o 

(2)n,JV=2 - (3)n,Af=2 
(2)n,JV=2 - (4)n,iV=2 




100 



500 



0.5 



-0.5 

-1 

-1.5 
-2 
-2.5 



(b) 



(3)n,jv=2 - (i)n,iv=o 
3)n,jv=2 - (2)n,iv=2 
3)n,jv=2 - (4)n,iv=2 
(3)n,jv=2 - (5)n,jv=4 

3)n,JV=2 - (6)n,AM 
(3)n,JV=2 - (7)n,AT=4 



100 

Atom-dimer distance R (a.u.) 



500 



Fig. 11. Same as Fig. II Ot but for two diabatic states of the 
77, iV = 2 manifold: (a) state (2), (b) state (3). 
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° 

N 

1 -0.5 

o 

z 

-1 





iii! i i 
(11)£ + ,JV=10 - (8)E + ,AT=8 
(11)Z + ,A'=10-(9)E + ,A'=8 

(ii)z + ,w=io- (io)s + ,w=io - 

(11)Z + ,W=10- (12)Z + ,Ar=12 
(11)L + ,JV=10-(13)r + ,AT=12 - 








, , , I . 


(a)' 



40 



100 



500 



(2)E",iV=2 - (1)E-,JV=0 
(2)E",AT=2 - (3)E",JV=4 




100 

Atom-dimer distance R (a.u.) 



500 



Fig. 12. Same as Fig. [JJJbut: (a) for state \ (11)S + ,N = 10), 
and (b) state \(2)£~,N = 2). 

trary, W pr increases with smaller 77. For given pairs of 
states, W pr even shows a resonant behavior, correspond- 
ing to a crossing of the diabatic curves. For example, 
we observe the crossings already mentioned above, be- 
tween |(3)27+,TV = 2) and |(4)27+, TV = 4), and between 
| (11)27+, TV = 10) and | (12)27+, TV = 12). We note the sharp- 
ness of the latter resonance compared to the former one, 
which confirms its more diabatic character. We also note 
two resonances between | (1)77, TV" = 0) and | (2)17, TV = 2), 
reflecting two intersections of those two curves at 7? = 
32 and 69 a.u.. Another surprising aspect is that of the 
|(3)77, TV = 2)-|(5)77, TV = 4) resonance, due to the fact that 
the numerator and the denominator of W pq (7?) both van- 
ish at 7? ss 50 a.u.. 

On each graph of Figs. [TUirT2l we also display two hor- 
izontal lines characterized by [W^rl = 10%, that we con- 
sider as the limit of applicability of the C n / R n expansion. 



This criterion gives different low- 77. limits 7?* according 
the state under consideration. It is quite similar for the 
states | (2)27+, TV = 2), |(3)27+,TV = 2) and | (2)77, TV = 2), 
respectively equal to 122, 100 and 102 a.u., and also quite 
close to the value R m — 102 a.u. determined in another 
way in Paper I. For state | (3)77, TV = 2), all the cross- 
ings take place in the same region, except one, with state 
|(4)77, TV = 2), for which 7?; = 261 a.u.. The very slow 
convergence of the corresponding W pr {R) function can 
be explained as follows. The two states, that belong to 
the TV = 2 manifold are characterized by a small energy 
difference: i.e. at 7? = 200 a.u., 3.3 x 10~ 4 cm -1 , com- 
pared to —5.4 x 10~ 2 cm -1 between (3) and (2), due to 
a compensation effect between the C5 and C§ coefficients. 
By contrast, they are still coupled at very long distance 
due to the dipolar interaction, equal to 4.4 x 10~ 4 cnr 1 
compared to —1.1 x 10~ 4 cnr 1 between (3) and (2) at 
200 a.u.. The particular occurrence can also be observed 
for higher TV, e.g. the couple | (6)77, TV = 4)-| (7)77, TV = 4), 
but it tends to damp as TV increases. Similar features have 
been noted in the case of two interacting atoms |47| . 



At last, Fig. rT2] illustrates two situations where the 
value 7?* is significantly smaller. On panel (a), we have 
plotted Wp r (R) for the state |(11)2T+, TV = 10). As shown 
on Fig. [31 the diabatic analysis is estimated to be valid 
down to 77* = 71 a.u.. The conclusion is similar for the 
|(2)27~, TV = 2) represented on panel (b), since 77* = 88 
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5 Conclusion and prospects terns and even to dimer-dimer systems, except hydrides, 



since the balance between the rotational and electrostatic 
energies does not strongly vary for other alkali- metals. For 
example, in the case of Li2+Li(2 2 P), the dimer rotational 
constant is about 60 times higher than that of Cs2, and 
the quadrupole moments and polarizabilities are smaller. 
But, since the electronic overlap becomes important at 
smaller distances, we can predict the existence of a cross- 
ing region, that we have estimated between 26 and 43 u.a. 



In this article, we have characterized the long-range in- 
teractions between an alkali-metal dimer in the electronic 
ground state and arbitrary vibrational and rotational lev- 
els, and an excited alkali-metal atom. More precisely, we 
have studied in details the interactions due to the perma- 
nent quadrupoles and to the induced dipoles of the two 
reactants. We have focused ourselves in a range of atom- 
dimer distances where the electrostatic energy competes 

with the rotational structure of the dimer, which allows 

i- , , ,, j-rr , , .. 11 i By contrast, the inclusion of the atomic fine structure 

couplings between the different rotational levels. J ' 

mi , , • i j i j , u ■ , will create a species-dependent situation. In the case of 

Ihe potential energy curves that we have obtained r r 

i, , . ,. c i j i i ii i_ lithium, the spin-orbit splitting of the 2 2 P state being 

are characteristic ot coupled-channel problems, showing m ' r r ° ° 

i .,j • n , , ,. T j 0.335 cm -1 , it is comparable to the typical rotational and 

particular avoided crossings. By calculating the Landau- ' 1 J * 

n uu-i-i-c rj-u ■ i. electrostatic energies, and may therefore complicate the 

Zener probability tor some oi those crossings, we have ° ' J * 

,, , ,, j , , , r j- , , . potential curves displayed in this article. On the contrary, 

seen that the dynamics can strongly turn from adiabatic v 1 J J ' 

,j.i,. i , i i ■ ,,. i , for the other alkali- metal atoms, the spin-orbit splitting 

to diabatic, when the dimer rotational quantum num- ' * 1 ° 

, . c , , , • ,. , , i , is so high (e.g. 17 cm -1 for Na(3 2 P)), that the two fine- 

ber increases. Such characteristics have already been pre- b 1 b v JJ ' 

dieted with the [OH] 2 complex ^SM and with polar l S structure components nP 1/2 and nP 3/2 will not be coupled 

i i ptti i -J.J. j j. i i i 4 c u by the electrostatic interaction. 

molecules [71] . submitted to an external electric held. J 

Because of the strongly attractive induced-dipole in- 
teraction, all the curves displayed in this article are at- 
tractive, which is an encouraging result in the prospect Acknowledgements 
of forming excited trimers by photoassociation. However, 
we note the absence of long-range potential wells, that 

can have a key role in the first observation of cold Cs2 Stimulating discussions with N. Bouloufa, V. Kokoouline, 
molecules |72j and R. Vexiau, are gratefully acknowledged. M. L. ac- 

We have illustrated our formalism in the particular knowledges the support of Triangle de la Physique in the 
case of a Cs2 dimer and a cesium atom. But we expect our framework of the contract QCCM-2008-007T Quantum 
conclusions to be generalizable to other atom-dimer sys- Control of Cold Molecules. 
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